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ABSTRACT 

With presently known input physics and computer simulations in ID, a self-consistent treatment 
of core collapse supernovae does not yet lead to successful explosions, while 2D models show some 
promise. Thus, there are strong indications that the delayed neutrino mechanism works combined 
with a multi-D convection treatment for unstable layers (possibly with the aid of rotation, magnetic 
fields and/or still existent uncertainties in neutrino opacities). On the other hand there is a need to 
provide correct nucleosynthesis abundances for the progressing field of galactic evolution and obser- 
vations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, 
i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treat- 
ment will not provide the correct yields because the effect of neutrino interactions is not included. We 
apply parameterized variations to the neutrino scattering cross sections in order to mimic in ID the 
possible increase of neutrino luminosities caused by uncertainties in proto-neutron star convection. 
Alternatively, parameterized variations are applied to the neutrino absorption cross sections on nu- 
cleons in the "gain region" to mimic the increase in neutrino energy deposition enabled by convective 
turnover. We find that both measures lead to similar results, causing explosions and a Yg > 0.5 in 
the innermost ejected layers, due to the combined effect of a short weak interaction time scale and 
a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak 
interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, 
and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these 
innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich 
environment results in enhanced abundances of *^Sc, "^^Ti, and ^*Zn as requested by chemical evolu- 
tion studies and observations of low metallicity stars as well as appreciable production of nuclei in the 
mass range up to A = 80. 
Subject headings: supernovae, nucleosynthesis 



1. INTRODUCTION 

The problem of core collapse supernova explosions is 
an old one and the attempts to understand the mech- 
anism have been ongoing for more than 30 years. The 
idea that a massive star proceeds through all burning 
stages from H to Si-burning, finally leading to the col- 
lapse of the resulting Fe-core to nuclear densities with 
formation of a ne utron star has long been discussed (e.g . 
Baade fc Zwickl Il 934: ODDcnhcimcr fc Snvdeil ITosa 
Arnett fc Schrammlil97.1|) . Since the 1960s the explo- 



sion mechanism has been relat ed to neutrino emission 
from the hot col lapsed core (e.g. iColgate k. Whit'3ll966l: 
iBethe fc Wilsonlll985t lBethelll99n|) interruptedW a pe- 
riod when it was speculated that the strength of the 
bounce at nuclear densities could permit shock waves 
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with sufficient energies t o lead to prompt explosions (e.g. 
iBaron et alJll985l Il987|) . However, introduction of pre- 
viously neglected neutrino scattering processes were in- 
troduced (e.g. neutrino-electron scattering), which per- 
mitted the replacement of lost low energy neutrinos, led 
to a continuous energy leakage and to th e death of the 
jprompt shock within 10 ms after bounce l|Bruennlll989t 
Mvra & Bludman 1989). 

Since then, and with the first neutrino d etection from a 
core collapse s upernova (SN1987A, see e.g. lKoshibalT99l 
lBurrowslll99ntl . the hope has been that further improve- 
ments would lead to successful explosions via energy de- 
position through neutrino and antineutrino captures on 
neutrons and protons (ve + n — > p-t-e~, S'e+P — > n-\-e^). 
Two different paths were explored. 1. Convective insta- 
bilities, but with still simplified neutrino transport, caus- 
ing either (a) convective transport in the core and lead- 
ing to higher neutrino luminosities fe.g. .Keil et al. 199^ 
or (b) highe r energy deposition efficiencies in convec- 
tive regions llMa.vle fc WilsonI 119881: iHerant et alJll994t 
IJanka fc Mulled 119961 iFrver k WarrenI I2004j) . 2. Im- 
proved neutrino transport schemes, leading to higher 
neutrino luminosities via the full solution of the 
Boltzmann transport e quation for neutrino scatterin; 
and neutrino react ions ijMezzacappa fc Bruenn>il993ai 
iMesser et al]ll998j) . 

However, in recent years ID spherically and 
2D rotationally symmetric radiation-hydro calcula- 
tions have not yet shown successful supernova ex- 



plosions wi th the present k n owled g e of physical pro- 
cesses. URamnn .Tankal 1200(11 IMezzar.anna, et 



cesses IIKa,mnn .lankal 1^1)1)11 IMezzar.anna 

200H iLiebendorfer et al.ll2001t iL iebendo rfer et alJ 



Buras et alJ 120031 iJanka et alJ 2003: Hix et all 
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Laneranke et alJl2003t iThompson et all 120031. A very 
recent simulation of a 11.2 M0 core collapse shows 
a chance for successful weak explosions i n a multi-D 
treatment with spectral neutrino transport ijJanka et al.l 
I2005|) . This leaves us with two dilemmata. First, the 
fundamental problem that the supernova mechanism is 
still not understood. Second, there seems no way to pre- 
dict the correct supernova nucleosynthesis yields. This 
is a problem in itself, but is also a limitation for the 
rapidly expanding field of galactic chemical evolution, 
which is being energized by the large amount of recent 
abundance observations from low metallicity stars (e.g. 



Argast et al.l200ll200llSneden et alEoollCavrel et all 
2004HHonda et al.l l2004>. 'Frebel et al."2005(l. 



20041: iHonda et al] l2004>. 'Frebel et al."2005 

Superno va nucle osynthesis prediction s have a long 
tradition (fWoosley & Weaver 1986; Thielemann^_et_^y 
19901 Imoslev fc Weaver. 1995; Thielemann et al.| 
199a iNomoto et al.l Il997t iN akamura ct al. 2001; 
Rauscher et a,1.l I2OO2I 



Chieffi fc Li^milgil I200I 
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Umeda fc Nomotol 



But all of these predictions relied on an artificially 
introdu ced explosion, either via a piston or a thermal 
bomb ijA^ifderheide et al"! I1991D introduced into the 
progenitor star model. The mass cut between the ejecta 
and the remnant does not emerge from the simulations, 
but has to be determined from additional conditions. 
While the usage of artificially introduced explosions is 
justifiable for the outer stellar layers, provided we know 
the correct explosion energy to be dumped into the shock 
front (on the order of lO^'^ erg seen in observations), 
it clearly is incorrect for the innermost ejected layers 
which should be directly related to the physical pro- 
cesses causing the explosion. This af fects the Fe-group 
composition, discussed in detail in |Thielem^^^^^ 
(^996), hereinafter TNH96, and iNakamurgTeTld 
l(l99 ^'l. which was also recognized as a clear problem 
by E hicffi & Limongi (2002) and Umeda & Nomoto 
(|2002). The problem is also linked to the so-called 
neutrino wind, emitted seconds after the supernova 
explosion, and considered as a possible source of the 
r-process to produce the heaviest eleme nts via neutron 
captures (Takahashi et al. 1994; Woosl ev et all 119941 
Qian et al. 1997; Thomoson ct al. 2001; Wanai o et al.l 
2001: Terasawa et al.. .2002: Thompson et al 2003il . 



An indispensable quantity to correctly describe nucle- 
osynthesis in the innermost ejecta is the electron frac- 
tion Ye = {Z/A) in the layers undergoing explosive Si- 
burning. This Ye is set by the weak interactions in 
the explosively burning layers, i.e. electron and positron 
captures, beta-decays, and neutrino or antineutrino cap- 
tures. The dominant reactions in hot photodisintegrated 
matter, consisting mainly of neutrons and protons, are 

+ P ^ n + . 

In section 0] we will show that these reactions lead to 
significant changes in Yf, toward equilibrium before the 
matter is ejected. We will also show that the resulting 
Ye in the innermost ejected layers is close to 0.5, in some 
areas even exceeding 0.5. This has been strongly postu- 



lated as a requirement in order not to violate abundance 
constraints from galactic evolution and solar abundances 
(TNH96). 

The question arises how one could realistically sim- 
ulate this behavior, given the existing problems with 
self-consistent explosions. Discussed improvements 
which could lead to successful supe rnova explosions 
are rotation an d magnetic fields (e.g. iThompsonI I2OOOI 
iThompson et ^1 . 2004) or uncerta inties in neutrino opac- 
ities (see e.g. iBurrows et al. ''20041 or other microphysics 
properties. They would introduce additional mixing at 
the neutrino sphere and convective transport or change 
the neutrino luminosity via improved opacities. This in- 
dicates the two options for successful explosion as dis- 
cussed above: (a) enhanced neutrino luminosities or (b) 
enhanced deposition efficiencies for neutrino capture in 
convective layers. These effects can be simulated in 
two ways: (a) Boosting the neutrino luminosity via a 
scahng (reduction) of the neutrino scattering cross sec- 
tions on nucleons while keeping the electron/positron and 
neutrino/ antineutrino capture cross sections on neutrons 
and protons at their original values, (b) Boosting the 
energy deposition efficiencies by enhancing the neutrino 
and antineutrino captures on neutrons and protons. Nei- 
ther approachrepresents a self-consistent treatment, but 
no external energy is required to produce (i) a successful 
explosion with (ii) a consistently emerging mass cut be- 
tween neutron star and ejecta. Moreover, our treatment 
guarantees that Ye is consistently determined by all weak 
interaction processes. 

2. MATTER IN A NEUTRINO FIELD 

Even though it is uncertain how significantly absorp- 
tions of the neutrinos emitted from the protoneutron star 
surface contribute to the revival of the shock, it is neces- 
sary to include the neutrinos and their copious interac- 
tion with the matter in the vicinity of the protoneutron 
star. If the explosion is launched such that the mass cut 
is directly determined by the early dynamics of the ex- 
plosion, this neutrino heated material will contribute to 
the deepest layers of the ejecta. If fallback occurs after 
the initial explosion, contributions by these innermost 
layers are only possible if strong mixing occures. In any 
case, this neutrino heated material will have significantly 
changed its composition with respect to its original pro- 
genitor composition. Hence, we investigate in this section 
the conditions in ejecta that are subject to large neutrino 
fluxes. 

The energy spectrum of the neutrinos is set in the 
vicinity of the neutrino spheres at a radius of initially 
~ 70 km. Before the launch of the explosion, about two 
thirds of the emitted neutrinos stem from the infalling 
matter which is squeezed in the gravitational potential 
and settles on the surface of the protoneutron star. The 
rise in the electron energies by the compression leads to 
more electron neutrino emissions than electron antineu- 
trino emissions. After the launch of the explosion, this 
contribution will decrease with the accretion rate and the 
less accretion-sensitive neutrino diffusion ffux from the 
hot protoneutron star will dominate. An energy equipar- 
tition among the different neutrino ffavors is expected to 
set in when the accretion luminosities have reduced to 
a negligible contribution. A more detailed description 
of this transition, however, requires multi-dimensional 
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simulations because the evolution of the accretion rate 
shows quite aspheric al features with na r row downflows 
and broad upflows fHerant et al.| [1994^ [ Burrows et aD 
[1995; Janka & Miillcr 1996: iBuras et alJl2003j) that are 
ignored in spherical symmetry. The emitted electron fla- 
vor neutrinos may essentially interact with the material 
behind the accretion shock out to radii of about 300 km 
( "essential" meaning electron fraction changes on a time 
scale of 100 ms). The interactions decrease steeply with 
increasing radius due to geometric l/r^ dilution of the 
neutrino field. 

In order to illustrate the basic behaviour of shock- 
heated matter in a neutrino bath, we first consider only 
the four dominant reactions, electron capture on free pro- 
tons e~ + p ^ n + i^e, and positron capture on free neu- 
trons e+ -|- n ^ p -|- and their inverse reactions. Two 
independent conditions are required to specify the elec- 
tron fraction and the entropy of the material, for exam- 
ple weak equilibrium and balance in the energy exchange 
with neutrinos. 

The change of the electron fracti on, Y ei with time, t, 
is given by Eqs. (C15) and (C20) in [Br]!Ien n (1985). The 
neutrino opacities, %, and emissivities, j, are linked by 
the reciprocity relation (detailed balance) described in 
Eqs. (C7) and (C8) in the above reference. The reci- 
procity relation involves the temperature, kT = 
the neutrino energy, E, measured in the rest frame of 
the fluid, and the chemical potentials, ^p, and fie, 
for neutrons, protons, and electrons respectively. We 
can therefore label contributions from electron, positron, 
neutrino, and antineutrino capture with EC, PC, NC, 
and AC respectively, and express the opacities in NC 
and AC by the neutrino emissivities. After having col- 
lected all terms that do not depend on the neutrino en- 
ergy into a common factor, K, we write the total change 
in the electron fraction in the following form: 

-'^=K ! dEE^ \h(E + Q) (-EC + NC) 
c at J 

+ e{E-Q-me) h{E - Q) {PC - AC)] . (1) 

Here, the details of the roughly quadratic energy depen- 
dence of the cross sections are hidden in the function 



h{x) 



m„c 



1/2 



and a step function 6 (a;) is used to describe the energy 
threshold in the positron and antineutrino capture reac- 
tions. A very similar equation can be used to describe 
the change of the specific internal energy, e, of the fluid 
due to neutrino interactions: 



1 de 
c dt 



= K j dEE^ [h{E + Q) {-EC + NC) 

+ e{E-Q-m„) h{E - Q) {-PC + AC)] . (2) 



The density in the supernova ejecta is low enough that 
we can neglect the nucleon degeneracy and nuclcon final 
state blocking described in Eq. (C14) of Brucnn (1985). 
The contributions to Eqs. (Q) and ^ from the individual 
reactions are given by 



NC = 

PC = 

AC^ 



1 

1 + g/3(-E-Q+A<e) 
g/3(£;-Q+Me) 



(3) 



where n„ and Up are the neutron and proton number den- 
sities, respectively, and f{E) is the neutrino distribution 
function in the rest frame of the fluid. 

Note that Eqs. (^ and ^ do not presume that the 
neutrinos are in equilibrium with matter, nor that they 
assume any particular spectrum. Some of these assump- 
tions, however, lead to useful analytical formulas for the 
equilibrium electron fraction. Dominance of the neutrino 
absorption terms has been assumed for the investigation 
of th e r-process in the neutrino wind of a protone utron 
star l|Qian fc Wo oslcv 1996: 'Thompson et alJ l200l1: and 
the cases where emission terms dominate or where the 
neutrinos are in thermal equilibrium have been analysed 
in a study of gamma-ray burst fireballs l)Beloborodovl 
[2003)- The balance between the four reactions in Eq. 
^ is determined by the ratio of the neutron and proton 
number densities and the exponential exp(/3[i?±(5=p/ie]). 
Its energy integral depends on a competition between the 
neutrino energy, E (in NC and AC populated according 
to the neutrino distribution functions /), the mass differ- 
ence between neutrons and protons, Q, and the electron 
chemical potential, /ig. Depending on the conditions, 
any one of these three quantities can assume a dominant 
influence on the balance in above reactions. 

For neutrinos at high energies compared to \Q — fif,\ 
the term exp(/3i5) is large and the neutrino absorption 
terms, NC and AC in Eq. dominate over neutrino 
emission. Hence, if the abundance of these high energy 
neutrinos is large, the equilibrium is determined by 
the balance between neutrino and antineutrino absorp- 
tion and therefore dependent on the unknown neutrino 
distribution functions. The rate of change of Ye is given 

by 

dYe 

~ ^u^n ^ Ye{\y^n + ^V^p)- (4) 



EC- 



1 



I ^ g/3(B+Q-pe) 



^P (1 - fu) 



Using Eqs. (64a) and (64b) of lQian fc WooslevI llT993^ 

for X^^n a-nd Xu^p and using the fact that neutrino and 
antineutrino luminosities are similar (upper right panel 
of Figure [5J it can be shown that Y^ > 0.5 is achieved 
pr ovided that 4 (to„ — m„) > (£77 — e^) (see also Fig. 5 
in iQian fc Woo slev (1996) and the discussion in §3 of 
'Hoff man et 310199611 

However, at earlier time, and as long as the eventual 
ejecta are close to the neutron star, all four reactions 
in Eq. ([3J) are active and the neutrino distribution func- 
tions are nontrivial functions of the accretion rate, the 
distance from the neutrino spheres, and the local weak in- 
teractions. Changes in the electron fraction are not only 
determined by the local neutrino fluxes and spectra, but 
also by the ability of the matter to accept captured elec- 
trons or neutrinos at the given conditions. Especially, 
when the electrons are degenerate the electron chemi- 
cal potential can dominate the exponential for average 
neutrino energies. In this case, exp {(3 {E -\- Q — /ig)) in 
Eq. ^ is small and exp {f] {E — Q + /ig)) becomes large. 
Hence, neutrino absorption, NC, and positron capture. 
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PC, are suppressed and the electron fraction decreases 
because of more prolific electron captures and antineu- 
trino absorptions. Balance is only established when the 
ratio between proton and neutron number densities has 
sufficiently decreased to compensate for the suppression 
of NC and PC introduced by the exponential. This leads 
to Up < Tin and an equilibrium electron fraction Yg < 0.5. 

Finally, in a plasma with nondegenerate electrons, 
the electron chemical potential becomes rather small so 
that the neutron to proton mass difference, Q, may ac- 
tually dominate the exponentials in the balance equa- 
tions, making exp (/3 (i? + Q — ^e)) in Eq. © larger 
and exp {f3 {E — Q + fi e)) smaller. Reversing the trend 
under degenerate conditions, NC and PC are favored 
and Up > Tin is assumed in equilibrium. For not 
too different neutrino and antineutrino fluxes and spec- 
tra, equilibrium will establish at > 0.5. Ac c ording 
to the analytical investigation in iBeloborodovl l|2003fl 
for the EC and PC reactions, this situation will oc- 
cur if the electron chemical potential fulfills the con- 
dition jjLe < Q/2. The larger binding energy favors 
protons over neutrons. High electron fractions below 
Yf, = 0.5 have been p redicted for supe rnova explosions 
ijFuller fc Meveilll995l:lThompsonll2000D . But recent su- 
pernova simulations with accurate neutrino transport 
have even exceeded the estimates, consistently finding 
values of Y^. > 0.5 in the vicinity of the mass cu t in cxplo- 
sion settings llLiebendorfer et al.l2001 :lBuras et al...200S: 
I .Tanka et al.1 1200,'^ iThompson et al.1 120041: iPruet et al.l 
120051 iFrohlich et alJl2005air 

3. HYDRODYNAMICAL SIMULATIONS 

The framework for this investigation are spherically 
symmetric simulations with implicit ge neral relativistic 
Boltzmann neutrino transport, sec Mczzc aopa fc Messed 
l)1999|) and iLiebendorfer et all (2004) for a detailed de- 
scription of the code agile- boltztran. It features 
a dynamically adaptive grid l)Liebendorfer et alJ 120021^ 
that concentrates grid points at the developing mass cut. 
The simulations are performed until the density drops to 
~ 10^ g/cm^ in the region of bifurcation between the 
ejecta and the remnant. At this time, the mass con- 
tained in radial mass zones is becoming very small and 
the run experiences ill-conditioned Jacobian matrices in 
the Newton-Raphson scheme. The simulations are then 
continued by an explicit hydrodynamic code until the 
temperature falls below T = 2x 10* K. This code employs 
an exp licit difference scheme similar to lColgate fc Whitd 
l)196fiD an d a simplified nucle ar reaction network as ex- 
plained in lBravo et a D (EMS). 

We use two different approaches to enforce explo- 
sions in otherwise non-explosive supernova models. We 
parametrize the neutral current neutrino scattering opac- 
ities on free nucleons with a factor ranging from 0.1 to 
0.7 and use a finite differencing^ that helps to artificially 
increase the diffusive fluxes in regions of very high matter 
density. The net result is a faster deleptonization of the 
protoneutron star such that the neutrino luminosities are 
boosted in the heating region. For the sake of compu- 
tational efficiency, this first series of parametrized runs 

^ According to Eg . (91) in ^Me zzacappa fc B ruennI [T993al) in- 
stead of Eq. (56) in ILiebendorfer et alJI2nn3r see sections 3.3.2 
and 4.1 in the latter reference for details. 



TABLE 1 

Name and properties of the discussed runs 



Run 


Parameter 


-Ecxpl [org] 


mcut [M0] 


ttj>o[s] 


^cnd [s 


A60 


60% scatt. 


0.24 X 10^^ 


1.585 


0.46 


0.64 


A40 40% scatt. 


0.78 X 10^1 


1.511 


0.27 


0.53 


A20 


20% scatt. 


1.24 X 10^1 


1.444 


0.20 


0.44 


BOS 


factor 5 


0.31 X 10^1 


1.586 


0.38 


0.53 


B07 


factor 7 


0.78 X 10^^ 


1.531 


0.26 


0.43 


BIO 


factor 10 


1.12 X 10^1 


1.509 


0.24 


0.40 



Note. — The parameter of series A specifies the percentage of 
neutral current interactions considered in the model. The param- 
eter of series B specifies the reduction of the heating time scale. 
The time after bounce where we had to stop the runs with neu- 
trino transport is displayed in the last column labelled by tend ■ The 
time of the first appearance of positive velocities is given in the col- 
umn tyyo. The mass cut mcut has been determined at the point 
where the total energy integrated from outside inwards reaches a 
maximum. The estimate for the explosion energy i?expl has been 
composed from the total energy of the unbound material between 
the masscut and the shock front at t^nd (mostly material that was 
in NSE) and a correction for the total energy of the bound layers 
ahead of the shock at progenitor composition. 



(series A) has been calculated with the lowest possible 
angular resolution, involving only inwards and outwards 
propagating neutrinos. However, all of these measures 
only affect the propagation of neutrinos in the model; 
the models are still closed and respect energy and lepton 
number conservation. We expect that series A represents 
a simplification of the phenomenology of supernovae 
that would be driven by higher neutrino luminosities 
than in the standard cases, for e xample different forms 
of protoneutron star convection llWilson fc Mavlelll993l: 
Keil et al.lll996t IMezzacaopa et alJ 119981: IBruenn et alJ 
2nnT or improvements in the uncertain nuclear matter 
physics. 

With progress in computer speed and code paral- 
lelization, we were able to perform simulations using 
standard resolution (6 angular bins, 12 energy groups) 
for the Boltzmann neutrino transport in the parame- 
ter study for series B. Series B also includes the weak 
magnetism corr ections in the neutrino cross sections 
ijHorowitd l2002|) . Explosions are enforced by multiply- 
ing the absorptivities and emissivities (i.e. the reaction 
rates for forward and backward reactions in + n ^ 
p + e~ and De + P <^ n + e~^) in the heating region 
by equal factors. This reduces the time scale for neu- 
trino heating without changing the important equilib- 
rium Ye and temperature. We hope to mimic with 
this approach a potentially increased heating efficiency 
in the heating region as it is expected in combina- 
tion with overturn in this convectively unstable domain 
llHerant et alll994l:lBurrows et al.ll995HJanka fc MtilleJ 
1996: ' Mezzacappa et al.lll998t iBuras et a\l\200^ . 

All models are based on a progenitor model with a 
ma in sequence mass of 20 Mq (Nonioto & Hashimot^ 
Il988l|) . The parameters in series A and B are chosen 
such that each series contributes with a barely exploding 
model, an extremely exploding model (in terms of param- 
eter range, the explosion energy itself seems to saturate 
around 1.2 x 10^"'^ erg), and a model with average param- 
eter setting. Important properties of the different runs 
are listed in Tabled 

Figure ^ presents an overview of the shock trajecto- 
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(a) (b) 




Mass Outside Masscut [M^^J Time [s] 

Fig. 1. — (a) Trajectories of the (accretion-) shock as a function of time. Scries A (sohd hnes) explodes because of artificially increased 
core diffusion luminosities. The neutral current neutrino opacity is used as parameter with 20% of the standard values (thin line), with 40% 
(medium line width), and with 60% (thick line). Scries B (dashed lines) explodes because of artificially accelerated neutrino absorption 
and emission processes in the heating region. The acceleration factors are 10 (thin line), 7 (medium line width), and 5 (thick line). Sec 
also the legend in panel (c) for an identification of the shock trajectories. The dotted line traces the trajectory of a representative fluid 
element located at 0.005 Mq outside of the estimated mass cut. Its electron fraction and specific energy can change by weak interactions 
on a time scale represented by the dot-dashed lines (reaction time scale as a function of radius). The lower branch belongs to the infall 
phase and the upper branch to the ejection. 

(b) Energy scales sampled at the representative fluid clement as a function of time. Thick lines show the fluid temperature (solid), 
the neutrino temperature (dashed), and the antineutrino temperature (dot-dashed). Thin lines represent the electron chemical potential 
(solid), the mass difference between neutron and proton (dashed), and the temperature that would be obtained after infinite exposure to 
the neutrino field (dotted). The electrons in the representative fluid element are degenerate during collapse and after shock compression. 
Neutrino heating and expansion during the ejection lifts the electron degeneracy and the neutron to proton mass difference becomes the 
dominant energy scale. 

(c) Electron fraction as function of mass when the runs axe stopped (0.4 — 0.6 s after bounce) with (A) reduction of the neutral current 
scattering opacities (solid lines) , or (B) enhancement of the reaction timescales (dashed lines) . The open circles represent the final electron 
fractions of run A40 after the freezout of charged current interactions. 

(d) Abundances sampled at the representative fluid element. Thick lines show the electron fraction (solid), neutrino fraction (dashed), and 
antineutrino fraction (dot-dashed). Thin lines represent electron fractions that would obtain after infinite exposure to the neutrino field. 
Only neutrino absorptions have been considered for the solid line and only emissions have been considered for the dashed line. The dotted 
line includes all charged current reactions. Panel (b) and (d) show that the electron fraction is first kept high by neutrino absorptions, 
later by neutrino emissions. The electron fraction at freezout is determined by competition between the neutrino interaction rates and the 
matter ejection timescale. 
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ries. Runs from series A are presented with solid lines 
and runs from series B with dashed lines. The legend in 
Figure^^ also applies to Figure^. In all runs, the accre- 
tion front stalls at about 100 ms after bounce at a radius 
between 180 and 300 km, depending on the parameters. 
The accretion front is slowly receding in the more opti- 
mistic models. Shortly before 200 ms after bounce, the 
accretion front moves outward again. There may still be 
some additional delay until the inwards drifting material 
behind the shock reverses its velocity and starts to ac- 
cumulate kinetic energy for the ejection. This happens 
at 199 ms after bounce for the fastest explosion (A20) 
and at 461 ms after bounce for the slowest run (A60). 
Bruenn's suggestion to locate the mass cut where the 
integrated total energy of all external material assumes 
a maximum agrees well with the actual bifurcation in 
the mass trajectories. The mass cuts, mcut, range from 
1.444 M0 to 1.585 M©. Reahstic 3D calculations where 
convection (responsible for the corrections applied to the 
weak rates in both series A and B) turned on in a delayed 
fashion could delay the explosions and lead to larger mass 
cuts. 

4. CONDITIONS OF MATTER IN THE VICINITY OF THE 
MASS CUT 

In the following, we trace a mass element in the exem- 
plary run A40. We choose a mass element that is 0.005 
Mq outside of the mass cut. The trajectory of this mass 
element is represented in Figure^ by a dotted line: At 
first, the element is falling into the gravitational poten- 
tial. After 200 ms it passes through the accretion shock 
at about 300 km radius and is instantaneously deceler- 
ated. A second phase of drifting around in the heating 
region follows until about 400 ms after bounce. Finally, 
the mass element is ejected to larger radii. 

Figure^ illustrates important energy scales along the 
trajectory. The dashed and dash-dotted thick lines at 
the top of the graph indicate the neutrino temperature 
for the electron neutrinos and antineutrinos respectively. 
They show a rising trend in the first half of the graph. 
This is because the protoneutron star shrinks and the 
neutrino spheres become hotter as they shift deeper into 
the gravitational well. The discontinuity at the crossing 
of the shock front stems from the Doppler shift when 
the mass element crosses the velocity jump at the accre- 
tion front. The change of the rise into a decline around 
t = 350 ms after bounce is due to the decrease of the 
accretion rate after the lauch of the explosion. Rising 
neutrino temperatures are resumed at a very small ac- 
cretion rate after t = 450 ms. 

With the full neutrino spectrum and abundances from 
the simulation and with the matter density as input, we 
calculate the equilibrium matter temperature along the 
trajectory according to Eq. (0) by requiring de/dt = 0. 
For consistency with the simulation, we have also in- 
cluded the charged current reaction s with nuclei a ccord- 
ing to the simple model described in lBruennI l)1985j) . The 
dotted line in Figure CJd shows the equilibrium tempera- 
ture of matter subject to the neutrino luminosities (the 
neutrinos themselves are not in thermal equilibrium with 
matter, their temperature is set in the vicinity of the neu- 
trino spheres where the matter temperature is higher). 
The lower part of Figure QJd shows the matter tempera- 
ture (thick solid line) and the electron chemical poten- 



tial (thin solid line). The electrons are degenerate in 
the cool infalling matter. The first little blip in the tra- 
jectory after t = 200 ms is due to the burning of the 
initial silicon layer to nuclear statistical equilibrium. It 
causes a slight rise in the temperature and decline in the 
electron chemical potential. The pronounced step up in 
both quantities is due to shock compression when the 
mass element hits the accretion front. During the drift 
in the heating region, we note a temperature increase 
towards temperature balance (dotted line) by neutrino 
heating. The onset of the explosion during this time 
also leads to an expansion and drop in matter density. 
Both effects work together to lift the electron degeneracy 
shortly before 300 ms after bounce (crossing of temper- 
ature and electron chemical potential lines). The evolu- 
tion during the third phase (ejection) is characterised by 
a density decrease. The weak interaction rates decrease 
and the temperature declines due to adiabatic expan- 
sion. The electrons stay nondegenerate and the electron 
chemical potential remains smaller than the neutron to 
proton mass difference (dashed thin line). In contrast to 
the electron-degenerate conditions found in past super- 
nova simulations that fail to explode, the expanding hot 
plasma under neutrino irradiation favors electron frac- 
tions that exceed 0.5 as discussed in section |5] 

The lower part in Figure^i shows the neutrino and an- 
tineutrino abundances with dashed and dash-dotted lines 
respectively. The variations are due to density changes 
rather than luminosity variations. The upper part of Fig- 
ure ^ shows the electron fraction from the simulation 
(thick solid line) and the equilibrium value determined 
by Eq. (dotted line). The dash-dotted line in Figure 
^ shows the reaction time scale as a function of radius. 
The upper branch belongs to infall, the lower branch to 
the ejection. Outside a radius of 600 km the reaction 
time scale is much larger than the dynamical time scale; 
during the drift phase of our mass element in the heat- 
ing region it assumes values around 50 ms. Thus, the low 
electron fraction during infall is mainly set by the progen- 
itor model. Before the shock front is crossed by the mass 
trajectory, the equilibrium Y,, is also low because many 
neutrons are bound in nuclei and not available as targets 
for antineutrino absorption. After the shock transition, 
matter is dissociated and higher electron fractions are 
favored. At first sight, the equilibrium electron fraction 
appears higher than expected at the given electron de- 
generacy. The reason are neutrino absorption rates that 
are by an order of magnitude larger than the neutrino 
emission rates at these moderate temperatures. The thin 
solid line shows the high electron fraction equilibrium as 
it would evolve if only neutrino absorption were consid- 
ered. The emission reactions alone favor a much lower 
equilibrium Y,, (thin dashed line) because there are only 
few positrons to capture under degenerate conditions. 
With the following rise of the temperature, however, the 
neutrino emission reactions (e.g. electron capture) gain 
weight with respect to the absorption reactions and the 
equilibrium Y^ correspondingly adjusts to lower values 
in the time window between t — 235 — 275 ms. But as 
the electron degeneracy is lifted with further tempera- 
ture increase and expansion, and the electron chemical 
potential dips below half the neutron to proton mass dif- 
ference, the emission rates start to join the absorption 
rates in favoring higher electron fractions (steep rise of 
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the thin dashed hne) . The equihbrium increases again. 
The descent at very late time is, as in the beginning, due 
to the reappearance of nuclei. The electron fraction in 
the simulation (thick solid line) can now easily be under- 
stood: At each time it evolves towards the equilibrium 
value for the combined reactions (dotted line) at the pace 
of the local reaction time scale. It freezes out when the 
mass element is ejected. Note that for an analytical esti- 
mate of the electron fraction in our application one would 
have to combine the approximations for neutrino absorp- 
tion rates in Eq. (64) in Qian fc Wooslev (,1996) with the 
a pproximation for neu trino emission rates in Eqs. (9-10) 
in iBeloborodovl |)200.1D and to consider the reaction time 
scale in order to find the correct freeze-out value in the 
Ye evolution. 

We find that all simulations that lead to an explosion 
by neutrino heating develop a proton-rich environment 
around the mass cut with Y^ > 0.5. This is illustrated 
by the electron fraction profiles shown in Figure Q];. The 
open circles denote the final (i.e. at T < 2 x 10^ K) 
electron fraction for the run A40. The mass scale is nor- 
malized to the respective mass cut. The different runs 
from series A show an almost identical electron fraction 
profile at the mass cut. The competition by the reaction 
and ejection time scale is not directly influenced by the 
different explosion parameters, i.e. the enhanced neu- 
trino diffusion at higher densities. The electron fraction 
profiles of series B, however, respond to the different re- 
action time scales set by the explosion parameters in the 
heating region. 

The electron fractions around and outside of rricut + 
0.1 Mq are still close to the progenitor values. Differ- 
ences in this region stem from the different locations of 
the mass cuts within the progenitor composition. It is 
important to note that the investigated region at the 
mass cut is highly unstable against convection because 
of a large negative entropy gradient. It is likely that 
the discrep ancies in Y^ are h eavily mixed on a dynamical 
time scale (jKifonidis et al.l l2003). We expect, however, 
that the Y^ remains high in an averaged sense (see also 
|Pruet et al. 2005). Moreover, matter blobs that leave 
the heating region in an environment of large convective 
turnover may still show qualitatively similar features in 
comparison with our spherically symmetric shells, be- 
cause the high electron fraction in the neutrino field is 
enabled by the discussed general features of expanding 
hot matter. We believe that the dependence on the de- 
tails of our different simulations is small. 

5. NUCLEOSYNTHESIS 

For the nucleosynthesis results presented here, we con- 
sider only the first few zones outside of the mass cut 
enclosing a few hundredths of a solar mass where values 
of Ye higher than 0.5 are achieved. To determine the 
final electron fraction Ye in supernova ejecta it is neces- 
sary to include neutrino absorption reactions on neutrons 
and on protons as well as electron and positron captures 
reactions. 

An example of the influence of the individual weak in- 
teraction contributions leading to Ye > 0.5 is given in 
Figure 13 (bottom right) for an exploratory study of one 
mass zone. Also shown are the neutrino luminosities L^, 
(top left) and energies Si, — (E^) / (E^) (top right). For 
this exemplary mass zone it can be seen in the upper right 



panel of Figure |2 that (ef — e,y) is always smaller than 
4(m„ — nip). There are several phases that can be iden- 
tified in Figure 121 and that have been discussed in section 
121 At early times (t < 0.3 s) matter is degenerate and 
electron capture dominates. At the same time matter is 
being heated by neutrino energy deposition and around 
i w 0.3 s the degeneracy is lifted (see upper panel of Fig- 
ure 12). At this time, the ratio between electron captures 
and positron captures significantly decreases and neu- 
trino absorption reactions start to dominate the change 
of Ye and, as discussed before, the average neutrino en- 
ergies favor Ye > 0.5. As the matter expands the density 
decreases, reducing the electron chemical potential. This 
results in positron captures dominating electron captures 
beginning around t w 0.3 s. From this time, the com- 
bined effect of positron capture and Ve absorption con- 
tributes to the final increase of Ye . 

Figure 121 shows the time evolution for a representative 
layer of models A40 (scattering cross sections on nucleons 
reduced by 40%) for the whole computational time. The 
final decline in the electron fraction Ye is due to /3-decays 
of the nucleosynthesis products. 

The position of the mass cut emerges consistently from 
the simulation as the region of bifurcation in which the 
density has dropped below ^ 10^ g/cm^. Based on the 
temperature-density profiles of all the matter in our hy- 
drodynamical simulation, the detailed nucleosynthesis is 
calculated in a postprocessing framework for the temper- 
ature range T > 2 x 10^ K. The extended nuclear reaction 
network used consists of 1072 nuclei with 1 < Z < 50, 
see Table 12 Neutral and charged particle reactions are 
taken from the recent R EACLIB compilatio n (contain- 
ing experimental rates bv lAngulo et al.Nl 999H C NACREl 
and t heoretical predictions bv lRauscherfc_Th ielemannl 
1)2000(1 ■ and are the same as used in lSchatz et al.. (,200L1 ). 
For the weak interaction rates ( electron/positron cap- 
tures and beta decays) the rates bv lFuller et al.l l(1982blaD 
are used for nuclei with ^4 < 45 (sd-shell). In the mass 
range 45 < ^ < 65 (jof-shell) the ext ended tabulation by 
ILanganke fc Martfnez-Pinedol 1(20011) is used. In Table 
O a detailed list of nuclei is given for which the above 
weak rates were utilized. The rates for neutrino and 
antineutrino captures on nuclei for the whole range of 
nucle i in the network are take n from a recent calcula- 
tion (Zinner fc Langankd 1200-^ . based on the random 
phase approximation calculations of tLanganke & Kolb^ 
l|200lL l2002'l. A complete list of nuclei for which neutrino 
and antineutrino capture reaction rates were included is 
shown in Tabled Hence, all weak interactions respon- 
sible for changes of Ye are taken into account in the re- 
action network, namely: neutrino/antincutrino capture 
on free neutrons and protons, neutrino/antineutrino cap- 
ture on nuclei, electron/positron capture, and 13^/(3^ de- 
cays. Neutrino scattering processes do no contribute to 
abundance changes and are thus not included in the re- 
action network used for postprocessing. Nevertheless, 
neutrino-induced spallation reactions can change the fi- 
nal abundances of some nuclei and we will investigate 
this in future calculations. 

Figure 01 shows the abundances after decay to stability 
of all nuclei for model A40 integrated over mass zones 
with Ye > 0.5, including in total ~0.011 Mq. For these 
mass zones, we are only concerned with the Fe-group 
composition. In Figure [SI integrated abundances after 
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Fig. 2. — Time evolution after core bounce of an ejected layer at 0.005 outside of the mass cut from a 20 Mq supernova progenitor. 
Top left: Luminosity of neutrinos and antineutrinos felt by a Lagrangian mass zone. Bottom right: Electron fraction Y^. Top left: Average 
neutrino energy eu = (E^) / (Ev) (thick solid and dashed lines). The difference in average neutrino energy, {e— — ev) (thick dot-dashed line) 
and four times the mass difference between neutron and protons (thin dashed line) are shown in the lower part. Bottom right: Individual 
weak interaction contributions leading to > 0.5. The individual contributions from neutrino/antineutrino captures and electron/positron 
captures are a factor ten larger than the total resulting dYe/dt. 



TABLE 2 
List of nuclei in reaction network 



Element 


Mass Range 


Element 


Mass Ran^ 


i;e Element Mass Range 


n 


1 


H 


1-3 


He 


3-6 


Li 


6-9 


Be 


7-12 


B 


8-14 


C 


9-18 


N 


11-21 


O 


13-22 


F 


15-26 


Ne 


17-34 


Na 


17-40 


Mg 


20-38 


Al 


21-40 


Si 


22-42 


P 


23-48 


S 


24-50 


CI 


26-42 


Ar 


27-44 


K 


29-48 


Ca 


30-50 


SC 


32-52 


Ti 


34-54 


V 


36-56 


Or 


38-58 


Mn 


40-62 


Fe 


42-64 


Co 


44-66 


Ni 


46-68 


Cu 


48-70 


Zn 


51-74 


Ga 


53-86 


Ge 


55-78 


As 


57-80 


Se 


59-84 


Br 


61-86 


Kr 


64-92 


Rb 


66-92 


Sr 


68-94 


Y 


70-96 


Zr 


72-98 


Nb 


74-100 


Mo 


77-102 


Tc 


79-104 


Ru 


81-108 


RH 


83-110 


Pd 


86-114 


Ag 


88-116 


Cd 


90-118 


In 


92-120 


Sn 


94-126 



Note. — Nuclear species used in nuclear reaction network for 
postprocessi ng. The mass range given indicates for each element 
the nuclei with the minimum and maximum neutron number. 



decay to stability are presented for model B07. In this 
model, the zones with Ye > 0.5 enclose '^0.0066 Mq. 
The positions of the mass cut and the explosion ener- 



gies for both models are given in Tabled The isotopic 
abundances (relative to solar values) result from postpro- 
cessing based on the hydrodynamical profiles and from 
employing the full nuclear reaction network including 
neutrino and antineutrino capture reactions. Note that 
unlike earlier supernova nucleosynthesis simulations, nu- 
clei beyond A = 64 are also produced in appreciable 
amounts, ranging in fact up to A = 80 or even be- 
yond, due to neutrino interactions with matter during 
the whole period of explosiv e proce ssing . This is dis- 
cussed in detail in lFrohlich et al.l lj2005hl) . These nuclei 
are mainly produced in the zones close to the mass cut 
where the electron fraction depends strongly on neutrino 
captures. For these mass zones relatively high entropies 
are attained: s/kb ^ 30-51. In nucleosynthesis terms 
this corresponds to complete explosive Si-burning with 
a strong alpha-rich freeze-out which also leaves a finite 
proton abundance (O.OOO7M0) due to Ye being larger 
than 0.5. The high proton abundance permits the on- 
set of an rp-process which, however, does not proceed 
too far in A as (due to the high entropies) the densities 
are too small. The abundances result from the accumu- 
lation of matter at the waiting-point nuclei ^''Ge, ^^Se, 
^^Kr, ''^Sr. After decay to stability they produce the 
high abundances of ^'^Zn, ^^Zn, '^^Qe, and ^^Se. A rel- 
atively high abundance of ^^Kr is also obtained. ^^Kr 
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TABLE 3 

Nuclei for which electron and positron captures are 
included 



Mass Number Nuclides 




Time [s] 



Fig. 3. — Evolution after core bounce for a representative ejected 
layer at 0.005 Mq outside of the mass cut for model A40. For this 
layer the time for cooling from Tg = 2 to Tg = 0.8 is about 6 s. 
Top: Temperature (solid line) and density (dashed line) evolution. 
Middle: Entropy s per baryon in units of the Boltzmann constant. 
Bottom: Electron fraction Ye. The final decline is due to /3-decays 
of the nucleosynthesis products. 




Mass Number A 



Fig. 4. — Abunda nces for model A40 relative to solar abun- 
dances lLodder6l2003l) . Two different calculations are shown: with 
neutrino-induced reactions in the network (filled circles) and with- 
out neutrino-induced reactions in the network (open circles). 



1 


n, H 


21 


Mg, Na, Ne, F, O 


22 


Mg, Na, Ne, F 


23 


Al, Mg, Na, Ne, F 


24 


Si, Al, Mg, Na, No 


25 


Si, Al, Mg, Na, No 


26 


Si, Al, Mg, Na 


27 


P, Si, Al, Mg, Na 


28 


S, P, Si, Mg, Na 


29 


S, P, Si, Al, Mg, Na 


30 


S, P, Si, Al 


31 


CI, S, P, Si, Al 


32 


Ar, CI, S, P, Si 


33 


Ar, CI, S, P, Si 


34 


Ar, CI, S, P, Si 


35 


K, Ar, CI, S, P 


36 


Ca, K, Ar, CI, S 


37 


Ca, K, Ar, CI, S 


38 


Ca, K, Ar, CI, S 


39 


Ca, K, Ar, CI 


40 


Ti, Sc, Ca, K, Ar, CI 


41 


Ti, Sc, Ca, K, Ar, CI 


42 


Ti, Sc, Ca, K, Ar 


43 


Ti, Sc, Ca, K, Ar, CI 


44 


V, Ti, Sc, Ca, K, Ar 


45 


Cr, V, Ti, Sc, Ca, K 


46 


Cr, V, Ti, Sc, Ca, K 


47 


Cr, V, Ti, Sc, Ca, K 


48 


Cr, V, Ti, Sc, Ca, K 


49 


Fe, Mn, Cr, V, Ti, Sc, Ca, K 


50 


Mn, Cr, V, Ti, Sc, Ca 


51 


Mn, Cr, V, Ti, Sc, Ca 


52 


Fe, Mn, Cr, V, Ti, Sc 


53 


Co, Fe, Mn, Cr, V, Ti 


54 


Co, Fe, Mn, Cr, V, Ti 


55 


Ni, Co, Fe, Mn, Cr, V, Ti 


56 


Ni, Co, Fe, Mn, Cr, V, Ti, Sc 


57 


Zn, Cu, Ni, Co, Fe, Mn, Cr, V 


58 


Cu, Ni, Co, Fe, Mn, Cr, V, Ti 


59 


Cu, Ni, Co, Fe, Mn, Cr, V 


60 


Zn, Cu, Ni, Co, Fe, Mn, Cr, V, Ti 


61 


Zn, Cu, Ni, Co, Fc 


62 


Ga, Zn, Cu, Ni, Co, Fe 


63 


Ga, Zn, Cu, Ni, Co, Fe 


64 


Ge, Ga, Zn, Cu, Ni, Co, Fe 


65 


Ge, Ga, Zn, Cu, Ni, Co 



Note. — All nuclei per given mass number for 

which weak interactions by Fuller et al. 1 1982b 3) a-nd by 
ILanganke fc Martinez-Pineda L2001) are used in the nuclear net- 
work. 



is considered to be produced by the p- or 7-process in 
the ONe layers of the st ar. Chemical evolution studies 
fe.g. lTimmes et ahlllQQSj) underproduce ^"^Zn by about a 
factor 5. A possible site for the production of ^''Zn is the 
modest early-time neutrino-driven wind occurring after 
core bounce in sup ernovae (Wooslcv & Hoffman 1992). 
lUmeda fc Nomotol l|2005i) have found that the '^'^Zn/^^Fe 
ratio is enhanced if Ye is close to 0.5 and the explosion 
energy is as high as ~ 10^^ erg. Our proton-rich en- 
vironment constitutes an alternative or complementary 
production site. 

For the intermediate mass elements the main improve- 
ment compared to earlier calculations is the higher pro- 
duction of individual nuclei like *^Sc and ''^Ti. Scandium 
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TABLE 4 

Nuclei for which neutrino and anti-neutrino captures are 
included 



Element Mass Range Mass Range 
capture) (77 capture) 



n 


1 


1 


H 


1 


1 


He 


6 


6 


Li 


7-9 


7-9 


Be 


8-12 


8-12 


B 


10-14 


10-14 


C 


11-18 


11-18 


N 


13-21 


13-21 


O 


14-22 


14-22 


F 


16-26 


16-26 


Nc 


18-29 


17-29 


Na 


20-32 


19-22 


Mg 


21-35 


21-35 


Al 


22-37 


22-37 


Si 


24-39 


23-39 


P 


26-42 


25-42 


S 


28-42 


27-42 


CI 


30-42 


29-42 


Ar 


32-44 


31-44 


K 


34-48 


33-48 


Ca 


36-50 


35-50 


Sc 


38-52 


37-52 


Ti 


40-54 


39-54 


V 


42-56 


41-56 


Cr 


44-58 


43-58 


Mn 


45-62 


45-62 


Fc 


48-64 


46-64 


Co 


50-66 


49-66 


Ni 


52-68 


51-68 


Cu 


54-70 


53-70 


Zn 


56-74 


55-74 


Ga 


58-78 


57-78 


Go 


60-78 


59-78 


Sc 


67-84 


66-84 


Br 


69-86 


68-86 


Kr 


71-92 


70-92 


Rb 


73-92 


72-92 


Sr 


77-94 


74-92 


Y 


79-96 


78-94 


Zr 


81-98 


80-96 


Nb 


83-100 


82-98 


Mo 


85-102 


84-100 


Tc 


87-104 


86-102 


Ru 


89-108 


88-104 


Rh 


91-110 


90-108 


Pd 


94-114 


92-110 


Ag 


96-116 


95-114 


Cd 


98-118 


97-116 


In 


100-120 


99-119 


Sn 




101-120 



Note. — Nuclides for which neutrino and anti-neutrino capture 
reactions are included in the nuclear network. The mass range 
given indicates for each element the nucleus with the lowest mass 
number and the nucleus with the highest mass number. 



is mainly produced by the /?^-decays originating from 
■*^Cr and "^^V decaying via ^^Ti to "^^ Sc. Different calcu- 
lations of abundance yie lds fTNH96. lWooslev fc Weaveil 
119951 iChieffi fc Limongi ..2002) fail to predict the ob- 
served abundance of scandium l|Gratton fc Snede'nlll991t 
iCavrel et al.ll200'^ . Our calculations show that Sc can 
be consistently produced with iron in the inner regions of 
the supernova where is higher than 0.5. The ejected 
yield of Sc is 10^® Mq which is a factor of 10 larger than 
the value obtained for a similar star by TNH96. If we 




Mass Number A 



Fig. 5. — Abundances for model B07 relative to solar abun- 
dances iLodders 2003). Two different calculations are shown: with 
neutrino-induced reactions in the network (filled circles) and with- 
out neutrino-induced reactions in the network (open circles). 



assume that our total production of Fe is similar to the 
one obtained in TNH96 our Sc yield will be consistent 
with observations. ^^Ti is underp roduced by a factor 5 
in the chemical evolution studies of lTimmes "etaP llTgg?). 
The nucleus ^^Ti originates from the decay chain of ''^Mn 
which decays via ''^Cr and ^^V to ^^Ti. After decay to 
stabihty, the resulting yield of ^^Ti is ~ 5 x IQ-^ Mq. 
We find that the origin of the differences in nucleosynthe- 
sis yields is a consequence of an electron fraction above 
0.5 which is due to a consistent treatment of all weak 
interaction processes on free nucleons. The obtained 
values are not sensitive to the inclusion of neutrino and 
antineutrino captures on nuclei. 

In the absence of a (yet) complete nucleosynthesis cal- 
culation covering the entire region responsible for Fe- 
group production we combine our abundances with the 
results of TNH96 (see Figures El and Ej) . The two cal- 
culations are combined in such a way that the resulting 
amount of Fe-group elements is the same as in this ear- 
lier work. For the inner zones, where neutrino and an- 
tineutrino capture reactions play an important role, the 
results of the present calculation are used. They consti- 
tute about 30% of the total production of Fe-group ele- 
ments. For the other zones where neutrino/antineutrino 
captures do not have significant influence on the final Y^ 
we use the abundance results from TNII96. This pro- 
cedure allows us to estimate the influence of a consis- 
tent treatment of weak interaction processes on the to- 
tal production of Fe-group elements. To further solidify 
these results full nucleosynthesis calculations are being 
performed based on the explodiri g mode ls. 

In a recent work, iPruet et al.l (l2005[) have studied a 
similar scenario (with a similar philosophy) for the in- 
nermost ejected layers, based on a 2D simulation by 
[janka et al. (2003). While the original 2D simulation did 
not yield an explosion, omitting the velocity-dependent 
terms from the neutrino momentum equation resulted in 
a successful explosion. They find results in agreement 
with our results. However, they do not report on the 
production of heavy nuclei with A > 6A due to the fact 
that they utilize the Y^ from the hydrodynamic calcu- 
lation but do not include neutrino interactions in their 
network. 
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Mass Number A 

Fig. 6. — Combined abundances of this work (model A40) and 
TNH96. The open circles are the combined abundances and the 
filled circles are the original abundances of the above reference. For 
details see text. 




55 60 65 70 
Mass Number A 

Fig. 7. — Combined abundances of this work (model B07) and 
TNH96. The open circles are the combined abundances and the 
filled circles are the original abundances of the above reference. For 
details see text. 



Figure |H1 shows elemental abundances of two calcula- 
tions compared to two s ets of observational data. One 
set of observational data ijGratton fc SnedenH igQll origi- 
nates from a sample of stars with —2.7 < [Fe/H] < —0.8, 
relevant for the average type II supernova contribution. 
The second set of observational data represents a s ample 
of extremely metal-poor stars ijCavrel et al.ll20d^ . The 
nucleosynthesis results of this work are shown in combi- 
nation with the results of the earlier calculation as shown 
in Figurc|51(thc results corresponding to Figurc[7|are very 
similar with the exception of Zn and heavier nuclei) . For 
comparison, the theoretical prediction by TNII96 is also 
shown. We see clearly an improvement for Sc and the 
heavy elements Cu and Zn. 

6. CONCLUSIONS 

Presently, self-consistent core collapse supernova sim- 
ulations in ID do not lead to successful explosions while 
2D models show some promise. Remaining uncertain- 
ties in neutrino opacities and/or the expected strong in- 
fluence of convection (due to hydrodynamic instabilities 
caused by entropy gradients and/or rotation and mag- 
netic fields) are likely to change this result. They may 




Ca Sc Ti 



Cr Mn Fe Co Ni Cu Zn 



Fig. 8. — Comparison of elemental overabundance in the mass 
range Ca to Zn for different calculations. The triangles with er- 
ror ba rs represent observat iona l data. The triangles facing up- 
wards IGratton fc Sneden^l^^9^^ originate from an analysis of stars 
with —2 .7 < [Fe/H] < —0.8. The traingles facing downwards 
ICavrel et al.i>20n3) is data for a sample of extremely metal poor 
stars (—4.1 < [Fe/H] < —2.7). The circles are abundances of this 
work combined with the work of TNH96 to obtain the same amount 
of Fe-group elements. The squares show the pure abundances of 
the previous reference. 



lead either to higher neutrino luminosities or higher ef- 
ficiencies of neutrino energy deposition via neutrino and 
antineutrino captures on nucleons. In order to exam- 
ine the accompanied nucleosynthesis in ID calculations 
of successful explosions, we performed simulations with 
variations in neutrino scattering cross sections on nucle- 
ons and/or neutrino and antineutrino captures on neu- 
trons and protons. In both cases successful explosions 
emerge with an interesting evolution of the Yg gradient 
in the innermost ejecta, which were followed up by a 
postprocessing for nucleosynthesis purposes. 

The detailed nucleosynthesis calculations with a consis- 
tent treatment of all weak interactions show an electron 
fraction Y^. > 0.5, i.e. a slightly proton-rich environment 
with relatively high entropies of up to ^50 ks per nu- 
cleon. This causes complete Si-burning with an alpha- 
rich (and proton-rich) freeze-out. About 0.0007 Mq of 
hydrogen remain in the innermost ejecta and do not stem 
from mixing this matter in from the hydrogen envelope. 
Such a proton-rich environment at relatively high en- 
tropies permits to produce also nuclei beyond A=64, up 
to A=80, with a major contribution to ^''Zn. The rp-like 
process does not extend to masses beyond A=80-100 as 
the high entropies imply too small densities for a path at 
very small proton separation energies. 

In addition, we find improvements within the Fe-group. 
The strong overabundances of ^*^^^Ni found in previous 
(too neutron-rich) environments are reduced. '^^Sc and 
^^Ti are enhanced to permit predictions closer to solar 
proportions. Especially the emergence of ^^Sc seems to 
be a solution to the previously not understood abundance 
of this only stable isotope of Sc. This discussion is also 
of interest with respect to ^^Ti, made in the alpha-rich 
freeze-out in the inner explosive ejecta. ^4^^ jg sensitive 
to Ye and reduced in the mass range where *^Sc is high. 
This will influence the overall predictions of '*''Ti. 

Values of Yg > 0.5 are due to the neutrino interactions 
with matter under electron non-degenerate conditions in 
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a convectively unstable domain and thus related to the 
explosion mechanism. The effect of neutrinos decreases 
with 1/r^ and about 50-60% of the Fe-group ejecta (the 
outer part of explosive, complete Si-burning) is deter- 
mined by values of Y,, equal to or close to the initial val- 
ues inherited from stellar evolution. In this first study 
we have tried to give an estimate for the overall Fe-group 
composition based on such a superposition of the present 
results for the innermost ejecta with those of TNH96. 
Future investigations will require to perform full nucle- 
osynthesis calculation for complete stars based on these 
exploding models. They will also require a sensitivity 
test of the nucleosynthesis results to the scaling factors 
for neutrino-induced reactions discussed in Figure ^ and 
in combination with the position of the mass cut and the 
explosion energy. This should be considered in order to 
reproduce results for supernovae where detailed observa- 
tional information in abundances, gamma-ray emitters 
and explosion energies is available. 
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